The shape of a flexible polymer in a cylindrical pore 
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Abstract 



We calculate the mean end-to-end distance (R) of a self-avoiding polymer encapsulated in an 
infinitely long cylinder with radius D. A self-consistent perturbation theory is used to calculate R 
as a function of D for impenetrable hard walls and soft walls. In both cases, R obeys the predicted 
scaling behavior in the limit of large and small D. The crossover from the three dimensional 
behavior (D — ► oo) to the fully stretched one dimensional case (D — ► 0) is non-monotonic. The 
minimum value of R is found at D ~ 0.46-Rf, where Rf is the Flory radius of R at D — > oo. The 
results for soft walls map onto the hard wall case with a larger cylinder radius. 

I. Introduction 

Beginning with the observation by Kuhn 1 that polymer coils are asymmetric even in dilute 
solutions, a number of studies have characterized the anisotropy of polymer chains 2 . Polymer 
molecules in good solvents, modeled using the Edwards model, are more anisotropic than Gaussian 
polymers because the number of ellipsoidal conformations in self-avoiding chains is far greater than 
spherical conformations^. The anisotropy of polymers, which is relevant in a number of applications 
involving polymer dynamics^, becomes even more pronounced in confined spaces. Nanopores 
(slits, cylinders and gels) align the polymer coils and distort their orientations, even when the 
characteristic confining volume is relatively large compared to the polymer volume*^. Confinement- 
induced alterations in the shape of a polymer is also relevant in biological applications. For example, 
a newly synthesized polypeptide chain transits the ribosome through a roughly cylindrical exit 
tunnel. The extension is perhaps achieved by an effective stretching forced f s ~ a ksT/D, where 
T is the temperature, ks is Boltzmann's constant, D (~ lnm) is the radius of the exit tunnel, 
and q is a constant. The magnitude of f s that is appropriate to the structure of the tunnel in the 
ribosome is between 4-10 pN, depending on a, which is large enough to unfold long stretches of 
proteins at low pulling speeds^. Another example is the encapsulation of a protein in the roughly 
cylindrical cavity of the E. Coli chaperonin GroEL-. In this case, substrate proteins are confined for 
a duration of time in a nanopore, which can enhance unfolding rates. In a very direct application, 
Tegenfeldt et. al.— have directly measured the genomic length of DNA molecules by trapping them 
in cylindrical nanochannels. 



Motivated in part by the above observations, we consider the behavior of a self-avoiding poly- 
mer of contour length L, confined to the interior of a cylinder of radius D. We are primarily 
interested in how the mean end-to-end distance, R, of the polymer changes as a function of D and 
the strength of the interaction between the cylinder and polymer. Daoud and de Genne&ii obtained, 
using scaling arguments, R when the interaction with the cylinder is purely repulsive. As D — > oo, 
the cylinder has no effect on the mean end-to-end distance R, which implies that R ~ Rf ~ IN U , 
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where I is approximately the size of one monomer, N is the number of monomers, and the Flory 
exponent v = 3/(cZ + 2) ~ 0.6 in d = 3 dimensions. As D — > 0, the polymer is effectively confined 
to d = 1. In the confined environment, there are only two relevant length scales, Rf and D, so 
that as D/Rp — > 0, and using the scaling assumption, we can write** 

R ~ R F f(R F /D). (1) 

As Rf/D — > oo, the chain is stretched in one dimension and becomes rod-like, thus resembling a 
one-dimensional self-avoiding walk. The scaling function f(x) takes the form 

where the unknown exponent m is determined from the condition R ~ iV as x — > oo, i.e. f(m+l) = 

I, so that m = 2/3 and thu&ii 

2 

R~i N (iy ( 3 ) 

The prefactor in eq. (j3J), which is a complicated function of and the polymer-cavity interac- 
tions, is difficult to compute. In this article, we calculate R for arbitrary values of D by adapt- 
ing the Edwards-Singh (ES) 12 uniform expansion method, which has been used in a number of 
application s W 3 . Note that, without the inclusion of an excluded volume term, the system will 
converge on a one-dimensional random walk, so that R ~ IN 1 / 2 , with no dependence on D. 

The remainder of the paper is organized as follows: In section II, we calculate R for a polymer 
confined to a cylinder with infinite polymer-cylinder repulsion. These calculations are repeated for 
soft walls in section III, and the differences between the two systems are determined. Finally, the 
effect of short ranged monomer-monomer interactions is briefly described in section IV. 

II. Hard Walls 

A self-avoiding chain is described by the Edwards Hamiltonian— /37~/[r(s)] =3/21 J Q L ds r 2 (s) + 
£>2, where B 2 = Vq Jq Jq ds ds' 5^(r(s) — r(s')) and Vq is the strength of the self-avoiding in- 
teraction. Following ES 1 - 2 , we replace the true Hamiltonian by a reference Gaussian, I3TLq = 
3/2/i Jq ds r 2 (s), where the effective step length Zi is determined by the ES method. We write 
(3H = PHq + £1 + B 2 where B x = 3(1// - l//i)/2 J Q L ds r 2 (s). We find^ =< r2 >= 
fdr dr L fV(r(s)) (r L - r ) 2 exp(-/3H[r(s)]) ~ < R 2 > -A + 0(Bf), with 

A = ^ < R 2 Bi > - < R 2 >o< Hi >o J + ( < R 2 B 2 > - < R 2 > < B 2 > ) , (4) 

and where < . . . >o denotes an average over the reference Hamiltonian /3Ho- The optimal value of 
l\ is chosen to satisfy < R 2 >=< R 2 >o, which is possible only if A = 0. The condition A = 
results in a complicated self-consistent equation for l\. For the unconfined case, ES showed 12 that 
higher order terms in the Bj's merely alter the numerical coefficient of l\ without affecting the 
Flory scaling laws. Thus, higher order terms will be ignored in this paper. For the remainder of 
the paper, all averages are taken with respect to the reference Hamiltonian 7Yo , so the subscripts 
on the brackets will be dropped. We also define Si =< R 2 > — < R 2 >< Bi >. 
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To calculate A, the Green's function for the reference Hamiltonian in a cylinder needs to be 
determined. Because of the infinitely repulsive polymer-cylinder interaction the Green's function 
vanishes at the walls of the cylinder. The resulting Green's function obeys the Heat Equation in 
an infinite cylinder, and the solution that satisfies the appropriate boundary conditions, in terms 
of the cylindrical coordinates r = (p, 0, z), is^ 



G(r ,r L ;L) = Gz{Z0 n2 ZL) £ f) cos m(fo - ^^uPo/D) J m (a p L /D) 

TT-LJ m=-oon=0 Jm+lV^mn) J m+l\Otm,n) 

where G z (zq,z£) = T>(z) exp(— 3/2Zi Jq ds z 2 ), J m (x) is the m th Bessel function, and a mn is 
its n th positive root. Using eq. < R 2 > becomes 



(5) 



< R 2 >= \Lh + ^ V (^-fl - 4-\- alJim ° 2 - -^e-< l ^ D2 \ ^ l -Ll l+ < R 2 >, 

(6) 

where N = ex P( — a o n hL /6D 2 ) / a^. The transverse term of the end-to-end distance, < R 2 . > 
=< x 2 > + < y 2 >, scales as < R 2 . >~ D 2 as D —> 0, implying that < z 2 >= i a L/3 ~ D 2 ^ as 
D -> (see eq. ©). 



Taking a derivative of ir givesJ* Si = — l/l\)d{R )/dl\, and we find 

1 5 r\i h)\ ^^\V < 2D 2 ) e 



(7) 



The second term in eq. (jlj is more complicated, but is simplified by splitting the averages 
into confined and unconfined terms. The unconfined z averages are calculated by completing the 
square in the exponent after Fourier transforming S(z(s') — z(s"))^^k. To compute S2 we define 

j 1 t i\ f 1 Jk( a kni x ) Jk(oCkn 3 %) J^rSS^mniX) , , oX 

Ik(m,{ni\)= / dx x— r- — — - — - — ^ -. r and (8) 



'0 Jk+l( a km) -4+l( a fcn 3 ) Jm+ 

E k (m,{ni};t,t') = I k (m,{rii}) exp ^ - ^^(a kni t < + a mn2 \t - t'\ + a fcm3 (l - i>))^, (9) 
where the 'time ordering' variables, i< and i>, are 

*<=<h , t>={ t ' do) 

with t = s'/L and t' = s" /L. In terms of these quantities, we find 



■Af V vr 3 /i ^ J Q J Q ^\t—p\\ a lni a lns 

' m,{rii} vii 

| £ (m,K};M') / _ _2 2_ _ < R| > fi£ ^ 

In the limit D — > 00, eq. converges to the Edwards-Singh self consistent equation for the 
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Figure 1: Plots of R/Rf as a function of D/Rp, for various values of L and Vo- The lowest curve 
corresponds to L/l = 1000 and Vq/1 = 0.1 and the uppermost curve has L/l = 1000 and Vq/1 = 0.5. 
There are two co-incident curves in the middle, one with L/l = 1000 and Vq/1 = 0.1, the other 
with L/l = 5000 and Vq/1 = 0.2. The minimum for all curves occur near D m i n ~ 0.46-Rf, with 

Rmin ~ 0.75-Rp. 

unconfined case^ 2 -, i.e. 



Thus, for large D, h ~ (24V r 2 / 2 L/vr 3 )i and R ~ R F ~ (2A/tt 3 )to (y Q l)§ L§. As L> -» 0, the ground 
state dominates in eq. (JSJ). Thus, only the (m, {nj}) = (0,0,0,0) = term in eq. (JHJ) makes an 
appreciable contribution. In this case, we find that 



For small cylinder radius, eq. JT3J) gives h ~ 0.393(y V^ 2 )^ an d ^ ~ 0.627(Vb//L> 2 )3L, after 

2 

Io(0) is evaluated. As predicted in eq. ©, i? has the proper scaling form ND~s. Thus, in both 
the D — > and D — > oo limits, the expected scaling form is recovered, including the predicted 
dependence of i? on D. 

To determine the behavior of R for intermediate values of the cylinder radius, the self-consistent 
equation (A = 0) can be numerically solved for l\ as a function of D, and R determined using 
R = (LZi)5. Fig. 1 shows R/Rf for different values of L and Vq. All of the plots are virtually 
identical, which implies that the only difference between the systems is the numerical value of Rf- 
As D/Rp exceeds unity, R — * Rf from below (fig. 1). In accordance with the scaling predictions, 
R/Rf increases substantially even under moderate squeezing (D/Rp ~ 0.2). Surprisingly, the 
crossover from the coil state (D/Rp 3> 1) to the stretched state {D/Rp <C 1) is non-monotonic. 
There is a minimum in R ~ 0.75Rf at D m i n ~ 0A6Rf (Fig. 1). This behavior, which has been 
previously observed for polymers confined to a slit^, is due to confinement-induced anisotropy 
in the polymer conformations. Monte Carlo simulations by van Vilet and ten Brinke 6 show that 
Rmin ~ 0.8Rp in a slit, which shows that confinement in a cylinder squeezes the polymer somewhat 
more than in a slit. 




(12) 




(13) 



III. Soft Walls 
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In real systems, the interaction between the polymer and the cylinder is not infinitely hard. 
It is therefore important to calculate R in the case of soft walls. By soft walls, we mean that the 
interaction between the polymer and the wall can be represented by a repulsive, non-hard sphere 
potential, which corresponds to heat conduction with radiative boundary conditions. In the hard 
wall case, G(r' , ro; L) = G(tl, r'; L) = 0, with r' = (D, <fr, z). If the walls are soft then the Green's 
function does not vanish at the boundaries, but satisfies^ dG s (r, r'; L)/dp\ p= n = CqG s {t, r'; L), 
where Co = oo corresponds to the hard wall case. The Green's function in this case is 

G?(t l ,t ;L) = -% £ Jk- cos ™(** ~ *>) Jm [ PO f™{ D) J ^lP™/D) e -p lnhL/ e D * ^ 

m __ oo n t> Iran Jm\Pmnj Jm\Pran) 

? (14) 

where we have defined the dimensionless parameters C = DCq, 7-mn = 1 + (Pmn ~ rn 2 )/C 2 , and 
where the /3 mn 's are the positive roots of 

When C> 1, (3 mn ~ a mn (l — 1/C), so that G s ~ (1 — 1/C)G by a Taylor expansion. It can then 
be shown that, for large C, R ~ (1 + 4/5C)i?. As D — > oo, the results for R for the hard and soft 
walls coincide. 



For finite C, the /3 mn 's can not be easily related to the a mn 's, so we define 

S " ^ < 2 >_ Afc ^ /3 2 n7o 2 n V 1 + C ~ (31 ) ' {3\ nl \ n K + C ) ' 

(16) 

The averages are easily computed by taking a derivative of eq. (|16f) . as in eq. (0). The £>2 
averages are tedious to calculate, but give 



, ^ (m,K};t,O ^ , 2 2 2 <R|> A) 

7o 2 ni 7o 2 n 3 V ^ ^ 3 2 ^ 6D2 ' 7/' ( } 

where Ef~(m, {rtj};i,i') is identical to Ey.(m, {rii}}t,t') in eq. @, except that a mn — ► (3 mn and 
I fc (m, {nj) -» I k (m, {nj), with 

Ik( m { n }) ~ [ dx x ^ k ^ kniX ^ Jk{flkn3 X ) Jm(l3mn2 X ) (^&) 
JO Jkifikm) Jk(Pkn 3 ) Jm(Prnn 2 ) 

As C — > 0, it can be shown that eq. (|15[) gives /?oo ~ v2C- In the small L> limit, for a finite 
Co, we find S\ converges to the first term in eq. f|13|). and 5*2 ~ — 8Vo/15Z) 2 (/iL/67r 3 ) 2 , implying 
Zi ~ 0.230(/Vb/-D 2 )sL and i? ~ 0.526(Vo?/-D 2 ) 3L. The scaling laws are unchanged by the softness 
of the walls; only the numerical coefficients are altered. 



The numerical solution for R as a function of D with varying Co shows that R m in ~ 0-75Rf 
as in the hard wall case (fig. 2). However, as Co decreases, D m i n decreases from ~ OAQRf to 
~ 0.32i?i? (Fig- 2). Thus, the soft wall behaves as a hard wall with a somewhat larger effective 
radius, D e ff = D + 5D. For D > D m i n , a shift of SD ~ D m i n (oo) — -D m , n (Co) causes R to coincide 
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Figure 2: Dependence of R/Rf as a function of D/Rp for various values of Co- The solid line 
corresponds to Cq/1 = 10 , the dotted line is for Cq/1 = 0.5, and the dashed line represents 
Co/l = 10 -5 . The values of L/l and Vq/1 are 1000 and 0.1, respectively 

with the end-to-end distance in the hard wall case. If we account for this shift, we find that the 
values of R for both cases differ at most by 5% for D > D m i n . Because the scaling laws change 
drastically for small D, a simple shift in D is not sufficient to reduce R to the value for hard walls 
for D < D min . 

IV. Effect of Monomer-Monomer Interactions on R 

Typically, there are interactions between monomers besides the universally present excluded 
volume interactions. As long as these interactions are short-ranged, the potential between mono- 
mers si and S2 can be modeled as B3 = —u)5(r(s2) — r( s i)) where As = S2 — si > 0. The 
insertion of this potential into eq. © leads to Vq — » Vq — uj5(s' — s{)S(s" — S2) = Vq — uS(t — 
ti)S(t' — t2)/L 2 in eq. (|11|). so that for infinitely hard walls, S3 ~ a;(3/87r 3 ZiAs)2 as D — > 00, 
and S3 ~ 2ujIq(0)/ D 2 (ZiAs/67r 3 ) 2 as D — > 0. Inclusion of monomer-monomer interactions for 
a polymer confined to soft walls yields the same scaling behavior in both limits, the only change 
being the numerical coefficients. Comparison of these scaling laws with eq. (|12|) and eq. (|13|) 
shows that there is an effective shift in the strength of repulsion due to a monomer attraction, with 
AVo ~ -wAs2i"2 for small D, and AVo ~ -wAs^^L^j for large D. In both regimes, the effect 
of 53 is insignificant compared to £2 for very long chains, provided that uj < Vq. Short-ranged 
interactions between all monomers can be computed by assuming that the potential is pairwise 
additive. In this case, given a distribution of interactions between the i th and j th monomers, u>ij, 
we find that Vq^Vq- J2 i;j u^<5 (3) (r(sj) - r(sj)), or V ->■ V - J dtdt' Vi(t, t')5(r(t) - r(t')), where 
Vi(s',s") = uj(s',s")/L 2 . Because S3 only produces a shift in Vq, the addition of this potential 
will simply reproduce the results in Fig. 1. for intermediate values of D. Thus, we expect the 
predicted scaling form to be unaltered for a heteropolymer confined to a cylinder, provided the 
effective intra-molecular interaction is not strong enough to induce chain collapse. 

V. Conclusions 

Inspired by a number of physical situations, we have calculated the dependence of the size of 
a polymer molecule confined to a cylinder using the Edwards-Singh uniform expansion method. 
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The theory presented here provides an approximate formula for the dependence of the end-to-end 
distance R for arbitrary values of the cylinder radius D. The major conclusions of the study are: 
(i) The theory yields, in the appropriate limits, the predicted scaling laws for R as a function of 
D. In particular, the expected scaling function is obtained in the D — ► limit. The advantage 
of the theory is that the numerical factors that are difficult to obtain using scaling arguments^ 
have been explicitly calculated. This allows for a calculation of R for any value of D. (ii) We 
have calculated R by numerically solving the self-consistent equation. Surprisingly, we found the 
crossover from the three dimensional behavior (D — > oo limit) and the fully stretched limit {D — > 
case) is non-monotonic. The minimum value of R ~ 0.75Rf is found at D ~ 0.46-Rf when the wall 
is infinitely hard. This is because the wall induces an orienting field that enhances the anisotropy 
of the polymer. This effect is greater for an encapsulated polymer in a cylinder compared to slit 
confinement, (iii) A direct calculation shows that asymptotic scaling laws are the same for both 
hard and soft walls. Any soft wall can be replaced by an equivalent hard wall with a larger cylinder 
radius, provided the wall-cylinder interaction remains short ranged. 
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